GRB Flares: UV/Optical flaring (Paper I) 
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ABSTRACT 

We present a new algorithm for the detection of flares in gamma-ray burst (GRB) light 
curves and use this algorithm to detect flares in the UV/optical. The algorithm makes use 
of the Bayesian Information Criterion (BIC) to analyze the residuals of the fitted light curve, 
removing all major features, and to determine the statistically best fit to the data by iteratively 
adding additional 'breaks' to the light curve. These additional breaks represent the individual 
components of the detected flares: Tgtart^ Tstop, and Tp^^ak- We present the detection of 119 
unique flaring periods detected by applying this algorithm to light curves taken from the Second 
Swift Ultraviolet/Optical Telescope (UVOT) GRB Afterglow Catalog. We analyzed 201 UVOT 
GRB light curves and found episodes of flaring in 68 of the light curves. For those light curves 
with flares, we find an average number of ^2 flares per GRB. Flaring is generally restricted to 
the first 1000 seconds of the afterglow, but can be observed and detected beyond 10^ seconds. 
More than 80% of the flares detected are short in duration with Ai/i of < 0.5. Flares were 
observed with flux ratios relative to the underlying light curve of between 0.04 to 55.42. Many 
of the strongest flares were also seen at greater than 1000 seconds after the burst. 



Subject fieadings: gamma-ray burst: general 

1. Introduction 

One of the many great advances made by the 
Swift (jGehrels et al.l 120041 ) mission was that of 
early time Gamma-ray Burst (GRB) afterglow 
follow up. Swift was specifically designed with 
rapid GRB afterglow observations in mind. Prior 
to the launch of Swift, GRB afterglow observa- 
tions generally did not start until hours after 
the burst, and an X-ray position was generally 
needed before any optical follow-up could occur. 
This meant that most optical detections did not 
take place until days after the GRB. Swift would 
solve this problem through the use of 3 separate 
instruments on the same spacecraft working to- 
gether. After the detection of a GRB by the 
Burst Alert Telescope (BAT; iBarthelmv et al 



20051 ) ■ Swift autonomously slews to the posi- 



tion, generally within ~ 100 seconds, allowing 



the X-ray Telescope (XRT; iBurrows et al]|2005bl) 



and UV/Optica l Telescope (UVOT;|Roniing_elal 
I2OOOU2OO4II2OO5I) to begin observations of the af- 
terglow. Swift has proven to be invaluable in fur- 
thering our understanding of GRB physics, having 
observed over 700 GRBs, and was specifically de- 
signed to be able to observe the early afterglow 
evolution and transition from the prompt emission 
to the afterglow stage. The early stages of the af- 
terglow proved to be very exciting and led to the 
discovery of a number of new featu res, including 



the " canonical" X-ray light curve (jNousek et al 
l200fi[ ). which has been observed in a number of 



GRBs (e.g. iHill et al.l I2OO6I lEvans et al.l l2009t ) 



-'Pennsylvania State Univ., 525 Davey Lab, University 
Park, PA 16802, USA 

■^Southwest Research Institute, 6220 Culebra Road, San 
Antonio, TX 78238, USA 

^MuUard Space Science Laboratory, University College 
London, Surrey, UK 



Another important feature seen ea rly in the Swift 
mission was X-ray flares (e.g. ([Burrows et al 



2005al: iRomano et al.|[2006l) ) 



Flares in X-ray light curves had been seen 
prior to their di s covery in XRT light curves (e.g. 



ery 

Piro et all Il998l \20m . but had only been ob- 
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served a handful of times. It was quickly shown 
that they are quite common, appea ring in approx- 



imate ly 50% of XRT afterglows (|0'Brien et al 



and are temporally displaced so as to be 
distinct from the prompt emission. Flares are ob- 
served as superimposed deviations from the un- 
derlying light curve and have been observed in all 
phases of the canonical X-ray light curve. 

Early in the Swift mission, several studies were 
performed that highlighted individual GRBs that 
exhibited either large numbers of flares or flares 
of unusually high fluence. Each of these stud- 
ies expanded our understanding of flares and 
the physical processes whereby they are cre- 
ated. In particular, t he studies of GRBs 050406 
(iRoinano et all l2006l). 050502B (iFalcone et al. 
20061), 050713 A (iMorris et al.l l2007l). 050724 



(ICampana et al. I2006F and 050904 (ICusumano et al 
20071 ) established the fact that flares are likely 
caused by internal shocks because of their steep 
rise and decay slopes, though the actual source 
of the flares is still debated and may be linked to 
instabilities in the ejecta or the release of stored 
electromagnetic energy. These studies also showed 
that X-ray flares are observed in both long and 
short GRBs, can contain energies as large as the 
prompt emission, they appear to come from a dis- 
tinct emission mechanism than the afterglow, are 
seen in both long and short GRBs, and can be 
temporally separated from the prompt phase by 
hundreds of seconds. Further studies have only 
reinforced these initial findings and have even 
shown that significant flares can be created at 
times greater than lO^s after the initial prompt 
detection (e.g. Swenson et al. 2010[ ). 

Some attempts have been made to look at larger 
collections of flares and have exami ned their prop- 



erties on a more generalize basis. iFalcone et al 



([20071) and IChincarini et al.l (j2007l ) examined the 
temporal and spectral properties, respectively, 
from a collection of flares found in 33 of the first 
110 GRBs observed by Swift. The combined re- 
sults from these two studies found that the late- 
time internal shocks were required to explain 10 
of the observed flares and that central engine ac- 
tivity was the preferr ed method for a r najor ity of 
the bursts. However. IChincarini et al.l (|2007() also 
state that more observations of flares over more 
energy ban ds are needed. A follow -up study was 
performed I Chincarini et al. 20091 ) that limited 



the data set to to those GRBs which had red- 
shifts, enabling a study of the actual energetics of 
the flares and found some indication that the flare 
energy may be correlated to the GRB prompt en- 
ergy, but was limited due to the number of bursts. 
They once again confirmed that more observa- 
tional work is needed. Additional studies further 
confirmed earlier results showing that X-ray flares 
are likely caused by late-time internal dissipation 
processes, which produces the prompt emission, 
and also showed that flares evolve over time, be- 
coming broader and flatter. However these studies 
limited their data to only the firs t 1000 seconds of 



the G RB afterglow light curve (jChincarini et al 



20101) or a li niited sample of 9 exce ptionally bright 



X-ray flares ijMargutti et al.ll2010t ). 

An attempt at incorporating inf ormation from 
multiple energy bands was made bv lMorrisi (|2008[ ) 
in which spectral energy distributions (SEDs) were 
created, using BAT, XRT and UVOT data, for 
flares f ound in the same s amp le of 110 GRBs 
used b y [Falcone et al. (2007) and Chincarini et al. 
(|2007h . The fits to the SEDs showed that the 



flares, unlike the afterglow, could not be fit by 
a simple absorbed power law. 

The number of studies analyzing flares in the 
UV / optical are even more limited than those for 



the X-ray (jRoming et al.l l2006a[ ) . The primary 



reason for this is the lower significance of most 
flares in the lower energy bands. While the X-ray 
flares are often easily identified by visual inspec- 
tion of the light curves, potential UV/optical flares 
arc more often overlooked or dismissed as noise. 

A notable example of flares detected by the 
UVOT i s found in the l i ght cur ve of the short GRB 
060313 iRoming et al.ll2006bl) in which late-time 



flaring was observed by the UVOT, but not seen 
in the XRT. The type of behavior seen in GRB 
060313 are not unusual to the study of flares in 
GRBs. The presence of flares in one bandpass, 
but not in another is seen often, and just as com- 
mon are cases where the flares are observed from 
X-ray all the way to the optical. In the specific 
case of GRB 060313 the flares could be consis- 
tent with density fluctuation in the circumstel- 
lar medium, provided that the cooling frequency, 
Vc, lies between the X-ray and UV/optical bands, 
which would explain why the flares appeared in 
the UV/optical but not in the X-ray. However, 
the flares can also be explained by central en- 
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gine activity at late times, which is the more 
preferred explanation for X-ray flares as stated 
above. Ai iother notable examp le is that of GRB 
090926A (jSwenson et alj 120101 ). This burst dis- 
played the previously mentioned late time flares 
at times greater than 10^ seconds, which can be 
explained by central engine activity at extremely 
late times, but also because the flaring is simul- 
taneously observed in the X-ray as well as the 
UV/optical. Identifying the source of the flares 
and whether X-ray and UV/optical flares have the 
same origin remains an important open question. 

The common factor in all of the aforementioned 
studies is that the flares were found by simple 
manual inspection of the light curves and were 
easily detectable by eye. This method has al- 
lowed for a significant number of X-ray flares to 
be detected, but has yielded a very small number 
of UV/optical flares due to the difficulties listed 
above. A blind, systematic search for flares in 
both X-ray and UV/optical bands has not yet been 
performed and is necessary to provide an unbiased 
sample of flares of all brightnesses. Such a sample 
would be able to address some of the limitations 
mentioned in the previous X-ray studies and would 
provide access to a relatively untapped source of 
knowledge with additional UV/optical flares. 

The complementary nature of two such flare 
catalogs would allow for more stringent constraints 
on the origin of flares in GRBs through cross- 
correlation of the two energy regimes. The pre- 
cise nature of the GRB central engine still remains 
largely unknown and, because flaring is most likely 
related to central engine activity, the study of 
flares is crucial to our unlocking of that mystery. 

In this paper we present the results from a 
blind, systematic search for flares in UVOT light 
curves. Using Monte Carlo simulations and dy- 
namic programming in conjunction with Bayesian 
statistics we have constructed the most complete 
catalog of UV/optical flares to date, and provide 
the temporal details of each flare, including Tp^ak , 
At/t, and the strength of the flare relative to the 
underlying light curve. In a forthcoming paper we 
will perform a similar analysis on the matching 
XRT dataset and present our joint analysis of the 
two catalogs. 

This paper is organized as follows: In §2 we de- 
scribe our data set as well as our methodology for 
identifying flares. We present our catalog of GRB 



flares observed by the UVOT in §3, and discuss 
the implications drawn from the catalog in §4. In 
§5 we present our future plans for utilizing this 
new flare detection algorithm. 

2. Methodology 

For the purposes of this study we will be us- 
i ng th e light curves produced by iRoming et al.l 
(|2013h . This Second Swift Uhraviolet/Optical 
Telescope GRB Afterglow Catalog provides a com- 
plete dataset of fitted UVOT light curves for both 
long and short GRBs observed by Swift from April 
2005 throug h Dec 2010, and mak es use of optimal 
CO- addition (jMorgan et al.l 120081 ) which results in 
better sampled light curves and increases the prob- 
ability of detecting flares. In addition to being 
optimally co-added, the Second UVOT GRB Cat- 
alog also normalizes the GRB light curves to a 
single filter from the 7 possible filters used dur- 
ing observations. This normalization further in- 
creases the timing resolution of the overall light 
curve during those periods when multiple filters 
were used during the same orbit. These two meth- 
ods combined, optimal co-addition and normal- 
ization, have resulted in a completely unique and 
previously unavailable data set that is suited for 
use in searching for fiares and other small features 
(as opposed to the First UVOT GRB Catalog; 
Roming et al. 20091 ). This sample is also approx- 



imately twice as larg e as the sample provided in 



Roming et al.l ( 20091 ). Due to the normalization 



of the light curves, we will not be performing any 
chromatic analysis on the light curves, and our ex- 
perience in fitting GRB light curves leads us to be- 
lieve that our detected fiares will not evolve signif- 
icantly with energy over the limited energy band 
observed by the UVOT. Our analysis described 
below is performed on the residuals of the fitted 
UVOT light curves, using t he fit ting parameters 
provided by Roming et al. ( 20131 ) when calculat- 
ing our residuals. 

Even with the increased probability of detecting 
flares that comes with using optimally co-added 
data, the flares that we hope to identify are below 
the significance level of the p revious X-ray stud- 
ies previously mentioii ed (e.g. iFalcone et al. 2007 : 
IChincarini et al. 2007) and we require a statisti- 
cally robust method to confirm that the fiares are 
real and not part of the background noise. For 
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this study wc have us ed the publicly av ailable R 
package struc change (jZeileis et al.ll2002[ ) and the 
brea kpoints function co ntained within the pack- 



age (jZeileis et al.l 120031) . The breakpoints func- 



tion specifically employs the approach of dynamic 
programming to compute the optimal number of 
breakpoints in a time series of data. In terms of 
GRB light curves, a breakpoint is defined as the 
time at which the subsequent data are no longer 
fit by the preceding slope (i.e. the beginning or 
ending of a flare, or any other break in the light 
curve). The computation of breakpoints involves 
breaking the time series into smaller pieces, itera- 
tively fitting the curve and minimizing the resid- 
ual sum of squares (RSS). As one would expect, 
simply minimizing the RSS would lead to fitting 
the light curve with n — 1 breakpoints (where n 
is the number of data points), with an individual 
segment connecting each data point. To counter 
this effect, the breakpoints function also c omputes 
the B ayesian information criterion (BIG; ISchwarz 
19781) ■ 



BIC = -2 X L + kln{n), 



(1) 



where k is the number of free parameters to be 
estimated, n is again the number of data points 
in the light curve and L is the maximized value of 
the log-likelihood function. 



L = ilogp{D\dj,M,) - logp(I?|^?,+i, A/,+i)), (2) 

which compares the probabilities of two possi- 
ble fits to the data, functions Mj and Mj+i, each 
with their respective set of parameters 0j and 9j+i, 
and returns the more likely fit based on the ob- 
served data points, D. The BIC is penalized and 
becomes increasingly large when either the data 
is poorly fit or when the number of free parame- 
ters increases and the data is overfit. The BIC is 
therefore minimized using the simplest model that 
sufficiently fits the data. The breakpoints function 
returns the optimal number of breakpoints by cal- 
culating the fit that minimizes the combined RSS 
and BIC values. 

The BIC is unlike many of the more commonly 
used statistical measurements used in astronomy 
(e.g. x^) F-test, etc.). The BIC is a "summary 
of the evidence provided by the data in favor of 
one scientific theory, represented by a statistical 



mode l, as opposed to another" (jKass fc Rafterv 
19951 ). but it does not provide a definite strength 
or probability t o a p referred statistical model. 
Kass fc Raftervl (|l995l ) do provide a guideline for 



interpreting the strength of evidence for the pre- 
ferred fit: 

BICi — BICmin Evidence Against Model i 
to 2 Not worth more than a 

bare mention 
2 to 6 Positive 
6 to 10 Strong 
>10 Very Strong 

For our purposes we will require BICi — 
BICmin > 6, or 'Strong' evidence, to determine 
the preferred fit. In the case of Bid — BICmin < 
6, we will adopt the simplest fit model (i.e. 
fewest breakpoints) that satisfies the criterion of 
BICi-i - BICi > 6. 

The breakpoints function does not take into ac- 
count the systematic and random error present 
in all data. In order to reintroduce the effects 
of these errors, we run a Monte Carlo simula- 
tion, randomizing the values of the observed data 
points in line with the measured errors. For each 
Monte Carlo iteration we calculate the optimal 
number of breakpoints and their corresponding 
times. For the purposes of this study we per- 
formed 10,000 Monte Carlo simulations on each 
GRB light curve. For each light curve we exam- 
ine the BIC value for each potential number of 
breakpoints (i.e. 1, 2, 3, 4...77, — 1) over all 10,000 
iterations and do the same for the RSS values. Ap- 
plying the BICi — BICmin > 6 criteria allows us 
to determine the statistically preferred fit to the 
data, that does not overfit the light curve. Each 
of the breakpoints found in this optimal fit corre- 
spond to a specific data point where the original 
fit to the GRB light curve no longer adequately 
describes the data and is likely caused by a flare. 
Appendix A provides a more detailed explanation, 
showing results of the breakpoints function run- 
ning on simulated light curves and flares. 

Once these potential flares have been identifled, 
we extract the properties of the flare from the light 
curve. A well deflned flare consists of three break- 
points: 1) Tstart, the time of the initial deviation 
from the underlying light curve decay, 2) Tp^ak, 
the time of the peak of the flare, where the slope 
goes from positive to negative again, and 3) Tstop, 
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the time when the decay of the flare returns to the 
underlying decay slope. In most cases we can not 
precisely identify the exact time of Tstart, Tpeak, 
or Tstop due to insufficient sampling of the light 
curve, or observing constraints creating gaps in 
the light curve. Because of these limitations we 
only calculate limits on the boundaries of Tstart 
and Tstop based on the available data, and cannot 
precisely define Tpeak, but rather use the observed 
data point exhibiting the highest flux during the 
period of flaring as a lower limit. Our estimates of 
the peak flux ratio, relative to the underlying light 
curve, will also be a lower limit due to the limi- 
tations in calculating Tpeak- This approach will 
ensure that we do not bias further studies with 
underestimated values of the peak flux. Our de- 
termination of At/t is also limited by the uncer- 
tainty in determining Tpeak because the true peak 
flux may occur any time between Tstart and Tstop, 
but because most flares are relatively short the er- 
ror in our estimate of At/t is only a few percent. 
In the few cases of flares with distinct features 
that could be analytically fit (e.g. Tstart and Tpeak 
for the fiare peaking near 80 kiloscconds in GRB 
090926 A shown in Figured]), these limits closely 
match the values derived from fitting the flare it- 
self. We therefore see no need to apply a different 
approach to the those few exceptional flares, by 
explicitly fitting them, but rather use the same 
limit approach as for the rest of the data set. 

Figure [T] shows the results of running the flare 
finding algorithm on the light curve of GRB 
090926A. Because the flaring occurs at such late 
times, the flares are stretched out and the individ- 
ual components become easy to see. The figure 
show the 9 identified breakpoints as vertical lines. 
These 9 breakpoints were then combined to cre- 
ate 3 individual flares, each with a Tstart, Tpeak 
and Tstop- The observing gap after the peak of 
the second flare means that we can only provide a 
lower limit on Tstop by placing it at the first data 
point after the gap. Using the fiare finding algo- 
rithm we were able to successfully i dentify the two 



flares that we previously identified (jSwenson et al 



20101 ). but also detected a third flare starting at 
the beginning of the light curve that we were un- 
sure of when attempting to identify flares by eye. 



3. Results 

Here we present the results of our analysis of 
the 201 UVOT GRB light curves from the Second 
UVOT GRB Gatalog. We find the presence of 
at least 119 unique fiares, detected in 68 different 
light curves. Table [T]) provides the following infor- 
mation for each flare: 1) GRB Name, 2) the flare 
peak time, defined as the observed time since the 
initial burst of the highest flux data point during 
the flaring period, as well as the limits on 3) Tstart 
and 4) Tstop, set to the last and first data points, 
respectively, that are well fit by the underlying 
light curve. 5) A limit on At/t based on the cal- 
culated peak time, Tstart and Tstop, and 6) a lower 
limit on the ratio of the peak fiux during the fiar- 
ing period, relative to the underlying light curve, 
using the observed fiux at the fiare peak time. The 
fiux ratio is normalized using the fiux of the un- 
derlying light curve to allow for direct comparison 
of each fiare across all light curves. Lastly, 7) the 
confidence measure of the detected flare as a frac- 
tion of the number of times the flare was identified 
during the 10, 000 Monte Carlo simulations. 

4. Discussion 

Our analysis shows that at least 33% of the light 
curves in the Second UVOT GRB Catalog contain 
possible episodes of flaring. This number is very 
much in line with analysis that has been performe d 
on X-ray light curves (e.g. Chincarini et al. 2010[) . 
This result, however, does not correct for those 
light curves that were so poorly sampled, or for 
which observations did not start until such late 
times as to make the detection of any flares chal- 
lenging. It is not unreasonable to assume that an 
even larger fraction of the light curves in the Sec- 
ond UVOT GRB Catalog contain flares that will 
simply never be detected due to these issues. 

For the purposes of this analysis we have di- 
vided the detected flaring period into three groups. 
First the 'gold' group, deflned as those flares 
with a confidence measure greater than 0.7 and 
a At/t < 0.5. This group provides a good detec- 
tion rate and satisfies the somewhat 'classic' def- 
inition of a flare. This group contains 46 flares. 
Next is the 'silver' group, which expands the pa- 
rameters to flares with confidence greater than 0.6 
and At/t < 1.0. This group includes the 46 'gold' 
fiares and adds an additional 24, for a total of 70 
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flares. The final group, the 'bronze' group, con- 
tains all 119 detected flares. 

For those light curves with detected flares, the 
most common number of flares for each group is 
1 per GRB with the average number of flares for 
each group being ^2. Fig[2]shows the distribution 
of flares per GRB for the gold, silver, and bronze 
groups, shown in black, blue and red, respectively. 
The flares in the gold and silver groups come pri- 
marily GRBs with a single flare, while the bronze 
group has a large fraction of its flares coming from 
GRBs with two or three flares. No GRB had more 
than three flares all belonging to the gold group. 
GRB 090618 displayed the most flaring activity 
with 6 unique flaring episodes detected, with all 
of its flares belonging to the bronze group. 

The earliest flare peak time occurs at 108 sec- 
onds after the trigger of GRB 060708, and the lat- 
est detected flare peaks at 787 kiloseconds after 
the trigger of GRB 050712, with 85% of all de- 
tected flares peaking before 1000 seconds. Figurc|3] 
shows the distribution of Tp^ak for the three groups 
of flares. The gold group has an average Tpeak of 
^^500 seconds, while the introduction of later flares 
in the silver and bronze groups push their aver- 
age Tpeak to ~900 seconds after the GRB trigger. 
However, the most commonly observed Tp^ak for 
both the silver and bronze groups is also ^500 
seconds. It appears that all three groups are sim- 
ilarly distributed, but with fewer late time detec- 
tions in the gold group due to the strict criteria for 
gold designation and the difficulties in detecting 
flares at late times. The distributions for all three 
groups can be reasonably fit by the same Gaus- 
sian, with different maximums, centered around 
~300 seconds with a full width at half max of 10 
seconds. 

The duration of the flares, taking into account 
the limited nature of our determination of T^tart 
and Tstop, vary from At/t of 0.01 to greater than 
10, with at least 80% of the bursts exhibiting a 
At/t < 0.5. Figure |4] shows the distribution of 
At/t for the three groups of bursts. As with Tpeak, 
all three groups can be flt by the same Gaussian 
function, with different maximums, centered at 
0.14 seconds and with a full width at half max 
varying between 0.17 seconds for the gold group 
to 0.21 seconds for the bronze group. The con- 
sistency between the three groups for both Tpeak 
and At/t shows that the detection algorithm is ro- 



bust even at lower confidence levels. It should be 
noted that for flares with At/t > 0.5, particularly 
those observed in the flrst few hundred seconds of 
the light curve, may act ually be the onse t of the 



forward shock emission ( Gates et al. 2009() 



The relative strengths of the flares varies from 
a minimum flux ratio of 0.04 to a maximum of 
55.42. Figure [5] shows the distribution of the flare 
flux ratios for the three groups. The bronze and 
silver groups can be fit by a Gaussian centered 
around ~0.5, whereas the gold would be better 
fit by a series of Gaussians due to the valleys be- 
tween the peaks at ^0.2, ^0.5 and ~1.3. This 
may indicate a preferred set of flare strengths in 
the UV/optical, though, a single Gaussian fit to 
the gold group does provide a centered value con- 
sistent with the bronze and silver groups, indicat- 
ing that it is likely a lack of data causing the poor 
flt and that all three groups are in fact equally 
well flt by the same Gaussian. A larger sample 
of gold flares, meaning continued GRB light curve 
observations in the UV/optical, will be required to 
determine whether it is a lack of data or an actual 
physical phenomenon causing the observed struc- 
ture. More than 83% of the bursts have flux ratios 
between 0.04 and 2, with representatives in each 
of the three groups, while there are 19 relatively 
strong flares with flux ratios > 2. Interestingly, 
14 of the 19 relatively strong flares are among the 
15% of flares that peak later than 1000 seconds. 
After the flrst 1000 seconds, light curves have gen- 
erally poor timing resolution due to the decaying 
nature of the afterglow, so it may be a simple ob- 
servational bias that leads to a large fraction of 
those late-time flares having large flux ratios (i.e. 
larger flares arc easier to to detect than small flares 
at late times). Further analysis and simulations 
will be required to determine whether an observa- 
tional bias exists or whether the high fraction of 
large late-time flares is indicative of a unique sub- 
set of GRBs capable of producing these types of 
flares. 

In addition to analyzing the burst parame- 
ters individually, we also performed an analysis of 
the UV/optical flare parameters compared to the 
GRB prompt parameters for each burst. Specifl- 
cally we compared the reported Tgo, prompt emis- 
sion fluence, and the amount of structure present 
in the prompt emission (i.e. single FRED-like 
peak versus multi-peak structure) to Tpeak, At/t, 
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the flux ratio and the number of flares per GRB 
and find no correlation between any of the prompt 
emission parameters and flare parameters. We in- 
terpret the lack of correlation to indicate that the 
emission source of the UV/optical flares detected 
is not the same as that of the high energy prompt 
GRB emission. 

5. Future Work 

We have introduced a new flare detection algo- 
rithm and have analyzed the GRB light curves in 
the forthcoming Second Swift Ultra violet /Optical 



Teles cope GRB Afterglow Catalog (jRoming et al 
20131) . We detect the at least 119 unique flar- 



ing episodes, many of which are previously unre- 
ported. The new method uses the Bayesian Infor- 
mation Criterion to determine the optimal number 
of breakpoints are needed to provide the best sta- 
tistical fit to the data. This method is generic in 
its application and can be used to find any sort 
of 'feature' which causes the light curve to deviate 
from the expected decay slope. We plan to apply 
this same detection algorithm to the matching set 
of GRB light curves as observed by the Swift X- 
ray Telescope and publish a catalog of X-ray flares. 
Using these two catalogs together we will correlate 
the X-ray and UV/optical flaring in an attempt to 
better understand the nature of GRB flaring and 
ultimately their origin. 
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A. Demonstration of breakpoints function on simulated data 

To verify the effectiveness of the breakpoints function, we performed a series of tests on simulated UVOT 
data. Starting with a simple power-law light curve, we induced a number of flares on the light curve and 
attempted to detect those flares using the method previously described. In order to test the ability of the 
breakpoints function to detect a wide variety of flares we varied the Tgtart time, the amplitude, the duration 
of the induced flares, as well as the slope of the underlying light curve. All of our analysis is done on the 
residuals of the best fit to the light curve. 



We use the criteria of BICi — BlCmin > 6, or 'Strong' evidence from iKass fc Raftervl (Il995l ). to deter- 
mine the appropriate number of breakpoints to assign to a specific light curve. We then attempt to group 
breakpoints together to form individual flares based on the relative position of the potential breakpoints 
relative to each other and relative to the underlying light curve. Our measure of confidence is determined by 
the number of times a specific breakpoint was identified in the course of the 10,000 Monte Carlo simulations, 
and also satisfied the criteria of BICi — BICmin > 6. 

As previously mentioned, the ideal case is when the amplitude and duration of the flare are both large 
enough to make T^tart , Tp^ak , and T^top easily identifiable. Figure |6] shows an example of simulated UVOT 
data with an induced flare that meets this criteria, and shows the potential breakpoints found by the code 
as vertical lines on top of the data. Analysis of the potential breakpoints in relation to each other and to the 
underlying light curve show that they collectively form a single fiare with the three potential breakpoints 
marking the approximate times of Tstart, Tpeak, and Tgtop- As with all the identified flares, the values 
assigned to Tstart and T^top are lower and upper limits of the actual values, with the last data point before 
the flare that is still well flt by the underlying light curve assigned to Tgtart and the flrst data point after 
the flare well fit by the underlying light curve assigned to Tstop- In the case of a well sampled light curve 
and flare, like that shown in Figure [51 the limits given by Tgtart and T^top do not differ significantly from 
calculated values derived from fitting the fiare. The difference is less by than the size of the exposure bin 
for data points we identify as Tstart and Tatop- We therefore see no need to differentiate our methodology 
and to fit those few flares that can be flt using a function. For this fiare, the measure of confidence for each 
individual breakpoint is very high. Specifically, the breakpoint associated with Tpeafc was identified in all 
10,000 iterations. In this case we are therefore 100% confident in the presence of a flare despite not detecting 
all three components at 100% confldence. We assign the overall confldence of the flare to be 10,000, reflecting 
this certainty. This example is truly the exceptional case, as UV/optical flares are rarely observed with such 
strength. 

A more typical size flare is shown in Figure [T] In this case, the vertical lines again show the positions 
of potential breakpoints found by the flare flnding code. Further analysis shows that these two breakpoints 
again form a single flare, however, only the approximate Tp^ak and Tstop times have been identified. A 
potential Tstart breakpoint was never identified during the 10,000 Monte Carlo simulations. This does not 
mean that the code failed to properly identify the fiare, but rather shows that even at low significance the 
code is able to detect fiares, however all the individual components of the flare may not be detected. In the 
case of the flare in Figure [3 the rise to the observed peak of the flare is short enough, due to a combination 
of a shorter flare duration and a small peak amplitude, that there are no observed data points during the 
rise of the flare. The first data point found by the flare flnding code that deviates from the underlying light 
curve is the observed peak of the fiare. In cases such as these the code picks the nearest observed data point 
before the observed flare peak as being the lower limit for the Tstart- In this case, the assigned time for 
Tstart is once again very close to the actual start of the flare. This flare is given a confldence measure of 
7,434, which is the confldence measure associated with the observed flare peak. 

In both of the previous examples, the sample of the light curves was continuous and uninte rrupted. How- 



ever, this is never the case with actual data. Our data from the Second UVOT GRB Catalog (jRoming et al 



20131 ) are influenced, at the very least, by the fact that the Swift satellite has a 96 minute orbit and that 
any target on the sky will be unobscrvablc for >50% of the orbit. Additionally any number of other factors 
including the observing of higher priority targets and conflicts with spacecraft observing constraints, have 
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produced light curves with uneven sampling and occasional large gaps. These gaps are the reason that we 
have decided to provide hmits on the values of Tstart and Tstop and to use the observed time of peak flux 
when reporting Tpeak- Figure [8] shows the resuhs from running the flare flnding code on the same basic hght 
curve and induced flares as in Figure [51 but with an observing gap overlapping the beginning of the flare. 
The code once again identifies three distinct data points that are potential breakpoints. Analysis of the data 
points shows that we are unable to determine whether they are all associated with a single flaring event 
due to gaps in the light curve. Not only are we uncertain whether the first and last potential breakpoints 
correspond to a Tstart and Tstop of a given flare or multiple flares, we are also no longer confldent that the 
point identified as the peak observed flux is an accurate approximation of the actual Tpeak or the peak flux 
level reached during the flare(s). Again we assign the last data point before the flaring period that was well 
flt by the underlying light curve as Tstart (in this case the last observed data point before the observing gap) , 
and do the same for Tstop by assigning the first data point after the fiaring period well fit by the underlying 
light curve (in this case the first data point after the second observing gap). Because we did not observe the 
majority of the flare and are uncertain in the actual peak flux achieved, the flux ratio we report becomes 
a lower limit, based on the peak observed flux and the flux of the underlying light curve at the time of the 
observed peak. This flare is given a confidence measure of 9,784, which is again the confidence measure 
associated with the observed peak. 

In each of these examples we have only shown cases where there are no additional potential breakpoints 
other than those associated with either T start i Tp^ak or Tstop- Occasionally the code does find more than 
three breakpoints for a single fiare, specifically in situations where there is poor data sampling during the 
peak of the fiare. In these cases, if no single data point stands out as a peak, the code will identify points 
to either side of the peak as being potential peak candidates, resulting in four breakpoints for a single fiare. 
In these cases we assign the data point with the larger fiux to be the peak flux time used in our calculations 
of At/t and the flux ratio. It should be noted, that this same series of four breakpoints could be the result 
of two individual flares occurring in quick succession. If the timing resolution is larger than the Ai/t of the 
flares, the code will not be able to correctly identify two individual flares, but will rather identify a single 
broad fiaring period. 
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Table 1 

Optical GRB flares from the Second UVOT GRB Catalog 



Source Name 


T) 


Tstart lower limit* 


Tstop upper limit* 


At/t 


Flux Ratio 


Confidence 








lower limit 




GRB050319 


333.58 


272.87 


428.21 


0.47 


0.83 


0.5075 


GRB050319 


1061.79 


927.06 


1208.41 


0.26 


1.19 


0.3259 


GRB050319 


802.84 


739.17 


927.06 


0.23 


0.31 


0.252 


GRB050505 


3181.79 


3181.79 


8065.87 


1.54 


2.11 


0.9408 


GRB050525 


271.40 


257.75 


299.17 


0.15 


0.56 


0.9644 


GRB050525 


608.92 


510.85 


636.67 


0.21 


1.24 


0.4149 


GRB050525 


186.92 


172.68 


214.72 


0.22 


1.03 


0.9651 


GRB050712 


787496.25 


714478.38 


844068.38 


0.16 


26.48 


0.9601 


GRB050721 


494.74 


423.97 


508.92 


0.17 


0.54 


0.4482 


GRB050801 


398.35 


383.12 


439.89 


0.14 


1.24 


0.885 


GRB050801 


341.22 


327.61 


355.35 


0.08 


0.72 


0.8037 


GRB050802 


887.98 


803.49 


1104.62 


0.34 


0.91 


0.5038 


GRB050802 


1428.47 


1299.30 


1492.69 


0.14 


0.70 


0.4371 


GRB050815 


131.78 


99.32 


145.33 


0.35 


1.64 


0.5691 


GRB050824 


94236.29 


81838.89 


203483.48 


1.29 


11.65 


0.7743 


GRB050908 


368.58 


214.89 


435.09 


0.60 


1.59 


0.9147 


GRB051117A 


610.93 


555.99 


773.26 


0.36 


0.85 


0.616 


GRB 060206 


129.90 


93.51 


1444.54 


10.40 


1.67 


0.6495 


GRB 060313 


524.12 


470.21 


795.27 


0.62 


1.39 


0.7843 


GRB060428B 


145182.38 


141352.78 


174712.64 


0.23 


8.90 


0.568 


GRB060512 


448.07 


432.01 


4019.50 


8.01 


2.25 


0.817 


GRB060526 


262.77 


242.76 


272.78 


0.11 


1.08 


0.9724 


GRB060526 


192.71 


177.44 


242.76 


0.34 


0.75 


0.9178 


GRB060526 


292.79 


282.79 


312.81 


0.10 


0.04 


0.5578 


GRB 060604 


203.71 


193.70 


213.57 


0.10 


0.44 


0.4402 


GRB 060708 


108.65 


98.65 


128.67 


0.28 


0.34 


0.6744 


GRB 060729 


18581.79 


12568.86 


31109.46 


1.00 


0.37 


0.9989 


GRB060904B 


256.86 


234.94 


284.31 


0.19 


2.13 


0.8475 


GRB 060912 


225.13 


209.74 


245.14 


0.16 


0.41 


0.7158 


GRB060912 


375.24 


365.23 


395.25 


0.08 


0.70 


0.5805 


GRB060912 


285.17 


265.16 


305.19 


0.14 


0.24 


0.5041 


GRB 061021 


222.52 


192.50 


232.53 


0.18 


0.54 


0.7412 


GRB061021 


4274.53 


452.69 


4683.61 


0.99 


0.31 


0.5813 


GRB061021 


5295.80 


4683.61 


5298.24 


0.12 


0.27 


0.5805 


GRB070208 


433659.09 


427659.38 


480022.66 


0.12 


1.52 


0.673 


GRB070318 


246.45 


226.44 


256.46 


0.12 


0.16 


0.2538 


GRB070318 


191458.19 


142816.11 


296319.50 


0.80 


2.63 


0.3125 


GRB070420 


1939.71 


1750.87 


5792.08 


2.08 


2.03 


0.6216 


GRB070518 


273.72 


243.70 


317.92 


0.27 


1.63 


0.8283 


GRB070518 


193.65 


177.69 


206.81 


0.15 


1.24 


0.5656 


GRB070518 


117.81 


107.80 


137.83 


0.25 


0.24 


0.4614 


GRB070611 


4733.36 


3347.02 


10492.06 


1.51 


0.78 


0.9438 


GRB070616 


1011.36 


846.23 


1149.26 


0.30 


0.90 


0.9267 


GRB070616 


787.52 


468.17 


816.63 


0.44 


0.47 


0.8939 


GRB070721B 


275.26 


255.24 


285.27 


0.11 


1.17 


0.9697 


GRB071031 


1118.38 


1102.01 


1158.87 


0.05 


1.18 


0.7647 


GRB071031 


576.25 


546.56 


842.05 


0.51 


1.32 


0.74 


GRB071031 


11666.85 


7596.43 


14239.37 


0.57 


1.60 


0.6842 


GRB071112C 


595.43 


572.52 


636.93 


0.11 


1.57 


0.7843 


GRB071112C 


245.42 


215.39 


285.45 


0.29 


1.38 


0.706 


GRB071112C 


18094.00 


12735.36 


45965.98 


1.84 


55.43 


0.6809 


GRB080212 


156.79 


126.10 


178.16 


0.33 


1.82 


0.7306 


GRB080212 


266.88 


223.14 


295.85 


0.27 


0.90 


0.535 


GRB080212 


356.95 


340.57 


378.72 


0.11 


1.09 


0.2891 


GRB080303 


573.73 


512.75 


620.90 


0.19 


2.07 


0.7797 


GRB080303 


223.47 


193.45 


256.43 


0.28 


0.76 


0.5921 


GRB080319B 


252.59 


232.58 


262.60 


0.12 


0.10 


0.6566 


GRB080319C 


1289.59 


1194.14 


1394.45 


0.16 


1.92 


0.8361 
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Table 1 — Continued 



Source Name 




Tstart lower limit* 


^sfop upper limit* 


At/t 


Flux Ratio 


Confidence 






(s) 






lower limit 




tjHrJUoUooU 


138.31 


128.30 


148.32 


0. 14 


0.46 


0.7335 




/I Oft OQ 


412.44 


453.43 


0. 10 


1.18 


0.712 


KjiXDVOUOZV 


192.45 


170.56 




0.84 


0.65 


0.4985 




146.84 


136.83 


166.86 


0.20 


0.26 


n Q o /I o 

u.oo4y 




123.54 


123.54 


133.55 


0.08 


0.08 


0.6278 


O D Ti>r\or\To 1 


300.46 


290.45 


330.48 


0.13 


0.26 


0.6111 




482.51 


412.46 


532.55 


0.25 


0.46 


0.8254 




1 13.06 


103.06 


133.09 


0.27 


0.16 


0.9201 


LiKJaUoUoiU 


229. 12 


199.09 


289.17 


0.39 


0.17 


0.7133 


K ±3 Uo Uy U O ±3 


286.85 


276.84 


306.56 


0. 10 


1.01 


0.5412 




507.00 


476.98 


527.02 


0.10 


0.36 


0.3713 


r~< T> 13 nsriQnf; 
LirtJjUoUyUD 


256.26 


241 . 19 


os/1 no 

Z(54.Uz 


0.17 


4.86 


u.oyoQ 


IjxlJjUWUy lo 


IZOoo.Di 


uUOZ. / z 


n /] r; QO o 
140(5y .ZO 


0.68 


2.49 


0. 7557 


D TU r\ o r\ m o 




OlZo4i .41 


OO ( 14o.ol 


0.68 


2.34 


0.5695 


LxKiaUoUy loA 


470.21 


460.21 


490.23 


0.06 


1 .26 


n "yom 


LixlJjUoUy iOA 


370. 14 


360.14 


390.16 


0.08 


1.47 


0.5295 


(jrtrSUoUyzo 


4ozy .UD 


zuyu. ( U 


AOAA Qn 

4y44.yu 


0.66 


1.18 


u.oyoo 


rlrS U o uy z 


247. 1 1 


217.09 


257. 12 


0. 16 


1.88 


0. 7826 




253.00 


one oc 

zUu.yo 


z / o.Oo 


0.26 


0.89 


0.6911 


OOTDnoi nnT 


433. 15 


413.13 


513.20 


0.23 


0.60 


0.5528 




nco no 


243.09 


302. 12 


0.23 


0. 14 


0. 7585 


r^T?Rnsi nriB 


iZDU.oD 


1 007 A 1 
iZZ ( .41 


1291 . 19 


0.05 


1.33 


0. 7367 


tjKiaUolUzy 


1 T1 o ct; 

lo / iz.bo 


CC C^O C Q 

ODOy .Do 


ooono A o 
zzUyU.4z 


0.85 


0.32 


0.8647 


/~i D T> no 1 lof: 


153.77 


143.76 


173.80 


0.20 


0.44 


0.7555 


tjKJaUyUizo 


1112.77 


950.67 


1 1 oo 

1 Izz.DD 


0.15 


0.53 


n oonT 
U.ooO ( 


^jrtJj Uy U i z tS 


ccQ no 
DOO. Uo 


608.41 


707.52 


0.15 


0.55 


0.8018 


LirtJjUyUizo 


1 -1 ns 
i4Uc5.0o 


1368. 14 


1467.41 


0.07 


0.42 


0. 7555 


xlJj Uy U4 U i Jj 


1117.76 


1 oco r^c 

lUbo.ob 


1187.56 


0. 11 


6.24 


0.8217 


LiKJaUyUoiU 


147.51 


132.47 


157.52 


0.17 


1.29 


0.7921 


O D 13 none on 


1 on A OK 


957.27 


1 c oo r^G 

looo.ob 


0.61 


2.97 


0.9025 


IjxlJj u y U o o U 


133.67 


123.66 


143.68 


0.15 


0.19 


0.8736 


rt Jj Uy U D o U 


173. 71 


153.69 


183.71 


0.17 


0.15 


0.8189 


OOTUnnncQn 
tjKoUyUDoU 


634.84 


595.07 


756.78 


0.25 


1.72 


0.689 


ODTUnnnci o 
CjKoUyUblo 


715.48 


705.59 


764.51 


0.08 


0.54 


0.8776 


OTDTUnnnci o 
tjKiaUyUDlo 


onon CO 

zusy .by 


oo"70 on 

zu /y .oo 


o 1 n /I in 
zlU4. lU 


0.01 


1.47 


0.7139 


ijrtJDUyUDlW 


ZOy / .4 ( 


n c Qo nn 


OfiOl QA 
ZO.il .04 


0.02 


0.95 


0.6723 


r^T?T3nonfii s 
Lirv.ljUyUDlo 


666.85 


641 .01 


691 . 17 


0.08 


0.68 


n anA ^ 


GRB090618 


1916.21 


1906.32 


1930.64 


0.01 


1.27 


0.55 


GRB090618 


829.04 


813.52 


853.39 


0.05 


0.44 


0.5474 


GRB090926A 


47722.57 


46791.30 


51749.93 


0.10 


2.51 


0.8781 


GRB090926A 


81701.3 


70136.27 


144308.98 


0.97 


1.12 


0.9078 


GRB090926A 


232412.66 


196569.48 


248216.91 


0.22 


0.65 


0.6123 


GRB091018 


341.45 


331.45 


371.48 


0.12 


0.26 


0.7297 


GRB091029 


1059.17 


1049.28 


1147.60 


0.09 


0.82 


0.6614 


GRB091029 


543.66 


523.65 


553.48 


0.05 


0.79 


0.5801 


GRB091127 


69405.64 


68961.88 


74007.65 


0.07 


1.80 


0.1644 


GRB100425A 


24126.33 


23640.73 


38746.47 


0.63 


2.82 


0.9208 


GRB100425A 


719.77 


651.76 


947.33 


0.41 


4.12 


0.8502 


GRB100805A 


151.24 


141.23 


171.25 


0.20 


0.35 


0.3225 


GRB100805A 


409.43 


392.57 


423.80 


0.08 


1.06 


0.3222 


GRB100805A 


630.03 


589.59 


707.20 


0.19 


1.56 


0.2733 


GRB100805A 


39387.52 


29555.01 


56840.63 


0.69 


1.97 


0.2561 


GRB100814A 


224.40 


214.39 


244.42 


0.13 


0.35 


0.3692 


GRB100814A 


284.46 


274.45 


294.47 


0.07 


0.46 


0.3079 


GRB100901A 


416.87 


396.86 


436.89 


0.10 


1.45 


0.8464 


GRB100901A 


228.79 


208.77 


238.80 


0.13 


1.55 


0.2395 


GRB100901A 


17679.91 


10070.17 


21634.66 


0.65 


0.27 


0.2295 


GRB100906A 


445.62 


424.48 


465.41 


0.09 


3.22 


0.7149 
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Table 1 — Continued 



Source Name 




T^tart lower limit* 


Tstop upper limit* 


At/t 


Flux Ratio 


Confidenee 




{^) 


{^) 






lower limit 




GRB101017A 


180.65 


144.90 


215.47 


0.39 


1.78 


0.7994 


GRB101117B 


164.51 


154.50 


174.52 


0.12 


0.47 


0.512 


GRB101117B 


307.99 


234.39 


328.00 


0.30 


0.62 


0.3753 



Note. — Flares are listed in chronological order by GRB date, then sorted by confidence. *A11 times are 
relative to the time of the initial burst trigger. At/i is calculated as {Tstop — T'start) /Tpeak- T^start and top 
arc lower and upper limits, respectively. Flux Ratio is calculated as the flux at the flare peak time divided 
by the extrapolated flux of the underlying light curve at the same time, normalized using the flux of the 
underlying light curve, and is a lower limit of the actual peak flux ratio. The confidence measure represents 
the fraction of times the flare was identified during the 10,000 Monte Carlo simulations. 



GRB 090926A 




10* 10^ 
Time Since Trigger (MET - 275631259) 



Fig. 1. — UVOT light curve of GRB 090926 A showing the 9 breakpoints detected by the flare finding 
algorithm. The 9 breakpoints are combined to form 3 individual flares. 
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Number of Flares Detected per Burst 
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Fig. 2. — Distribution of the number of detected flares per GRB. The three distributions are the gold (black 
font), silver (blue font) and bronze (red font) distributions described in the text. 
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Distribution of T 
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Fig. 3. — Histogram of the distribution of Tp^ak- The three distributions are the gold (black font), silver 
(blue font) and bronze (red font) distributions described in the text. Note: In cases where no additional 
flares added when going from the gold to silver group (or from silver to bronze) the bin of the former group 
is shown as a thin line on top of the latter group. 
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Distribution of A t/t 
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Fig. 4. — Distribution of At/t for the detected flares. The three flares with /S.t/t > 2.0 are omitted for 
scaling purposes. The three distributions are the gold (black font), silver (blue font) and bronze (red font) 
distributions described in the text. Note: In cases where no additional flares added when going from the 
gold to silver group (or from silver to bronze) the bin of the former group is shown as a thin line on top of 
the latter group. 
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Distribution of Flare Flux Ratio 




2.0 2.5 3.0 
Flux Ratio 



5.0 



Fig. 5. — Distribution of fiarc flux ratio, relative to the underlying light curve. The flares with flux ratios 
> 5 are omitted for scaling purposes. The three distributions are the gold (black font), silver (blue font) and 
bronze (red font) distributions described in the text. Note: In cases where no additional flares added when 
going from the gold to silver group (or from silver to bronze) the bin of the former group is shown as a thin 
line on top of the latter group. 
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Simulated UVOT Data 
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Fig. 6. — Simulated UVOT data showing an optimal flare. All three components of the flare, Tstart, Tj 
and Tstop, are detected and identified. 
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Simulated UVOT Data 
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Fig. 7. — Simulated UVOT data showing a flare with a small amplitude and duration. Due to the abrupt 
rise to the peak of the flare, the code does not identify a unique point as being associated with the start of 
the flare. We assign the first point prior to Tpeak to be Tstart- The data point assigned as Tstart is identified 
by the red line. 
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Simulated UVOT Data 
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Fig. 8. — Simulated UVOT data showing the influence of observing gaps on the flare finding code. An 
elevated flux level is identified as a potential Tp^ak, however the beginning and ending of the flare are not 
observed. The first data point prior to the observing gap is designated as T^tart and the first data point 
after the second observing gap is designated as Tstop- 
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